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^- ■ Anisotropic London theory is used to predict the Abrikosov flux lattice 

for arbitrary field orientations in crystals of YE^C^C^.-^ by minimising the 

Gibbs free energy for samples of a planar geometry. At low fields the "vortex- 

^" 

ON . chain" state exists, i.e. the inter-chain vortex separation scales as the inverse 



flux density, 1/B, whereas the intra-chain vortex separation is constant. At 



i-rt ■ higher fields there is a cross-over to the uniaxially distorted hexagonal lattice 

where the inter- vortex separation scales as 1/B 1 ' 2 . At low fields the vortex 



lattice is inclined towards the c axis, rotating towards the applied field at 
higher fields. The results of this calculation are in close agreement to the 
Bitter pattern experiments of Gammel et alB assumming an anisotropy ratio 
7 of 5 and an in-plane penetration depth A a of 1 413 A, indicating that the 
"vortex-chain" state has been observed in these experiments. 
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I. INTRODUCTION 

An interesting aspect of high temperature superconductors is the nature of the mixed 
state in which magnetic flux lines (or vortices) penetrate the superconducting sample when 
the applied field exceeds H c i. In isotropic materials the vortices are arranged in a regu- 
lar hexagonal latticeS. In the anisotropic high temperature superconductors, however, the 
vortex lattice is distorted for an applied field not parallel to the crystalline axis. The de- 
tailed structure and arrangement of vortices in these materials has been the subject of much 
experimental and theoretical attention in the last few years. 

For high magnetic fields the distortion of the equilibrium flux line lattice depends on the 
anisotropy ratio and orientation of the lattice with respect to the crystalline axisH. At low 
magnetic fields, however, it has been predicted that the flux line lattice becomes distorted 
into a series of equally spaced vortex chains, running parallel to the plane containing the 
applied field and the c axis: the 9 plane. This "vortex-chain" structure arises from the 
vortex-vortex interaction becoming attractive for magnetic fields not parallel to one of the 
principle crystalographic axesaJl The minima of the interaction are situated at a distance 
x m i n on either side of the vortex core in the 9 plane. The lattice is no longer expected to 
scale uniformly with field; the intra-chain spacing, the distance between vortices within a 
chain, is determined by x m i n . For a fixed magnetic flux density the intra- vortex spacing re- 
mains constant, and hence the inter-vortex spacing, the distance between chains, is inversely 
proportional to the flux density. 

A number of experiments have investigated the flux lattice directly; both on the surface 
of the sample, and within the bulk. Small angle neutron scattering experiments have been 
performed to probe the magnetic field distribution within single crystals at high fields. 
Yethriaj et alB observed a distorted hexagonal lattice in YE^CusOr-^, in agreement with 
the predictions of London theory. Keimer et alB, however, observed vortex-chains for an 
applied field of B = 0.5T inclined at 80° to the c axis. The presence of vortex-chains at these 
large fields is attributed to an exponential softening of the vortex lattice shear modulus, as 



predicted by Ivlev and KopninO. 

High-resolution Bitter decoration experiments have also been performed to directly 
observe the flux line lattice emerging from the surface of single, twin free crystals of 
YBa2Cu307_<5 in low fields by Gammel and et aW. Their results show the variation of the 
intra- and inter- chain distance as a function of the normal component of the applied field, 
Hz- These results indicate that both the intra- and inter- chain distances scale as roughly 
1/yH Z-, and hence do not appear to obey the vortex-chain state predictions. However, 
since the experiments are performed on flat, platelet crystals, demagnetisation effects play 
an important role in determining the density and orientation of the flux lattice with respect 
to the crystalline axis. It is the purpose of this paper to show that once the demagnetisation 
effects have been considered the experimental results fit the vortex-chain predictions very 
well. This implies that, although the Bitter decoration experiments are surface probes, they 
do provide useful information on the bulk properties of the lattice. It also demonstrates that 
three dimensional anisotropic London theory provides a good description of the flux lattice 
in YBa 2 Cu307_,5. 

The plan of this paper is as follows: in §11 we discuss the thermodynamics of planar 
superconductors in a magnetic field. This is the situation relevant to most of the recent 
decoration experiments on high temperature superconductors. §111 discusses the magnetic 
field distributions obtained from anisotropic London theory, which enables the pairwise 
interaction between vortices to be calculated. In §IV we explain the details of the lattice 
summation, for both the low and high field limits, while §V describes the results and a 
comparison with the experimental observations. We summarise in §VI. 

II. THE THERMODYNAMICS OF PLANAR SUPERCONDUCTORS 

The high resolution Bitter pattern experiments of Gammel et alBwere performed on flat 
platelets approximately 0.5mm x 0.5mm square, and between 5 to 40 /zm in thickness. Thus, 
demagnetisation effects must be considered in the determination of the flux line lattice. In 



figure 1 we illustrate the coordinate system used in this paper. The (X, Y, Z) axes define 
the crystallographic axes with Z being parallel to c and (X — Y) lying in the (a — b) plane. 
In general an applied field, H a , will be oriented at an angle <fi with respect to the c axis. 
Since we are only concerned with uniaxial anisotropy, the plane enclosed by c and H a is a 
symmetry plane. Hence, the magnetic flux density, B, defined as the average flux density 
per unit cell, also lies in this symmetry plane, at an angle 9 with respect to the c axis. We 
therefore define a coordinate axes (x, y, z) by a rotation of 9 about the Y axis, such that z 
is parallel to B. 

To simplify matters we will assume that the platelets are infinite in the a — b (X — Y) 
plane. The boundary conditions are then that the normal component of B and the parallel 
component of H to the surfaces are continuous, i.e., 

B z = B a cos<p = Bcos9 (1) 

and 

# x = # a sin0, (2) 

where B a = H a . Equation (1) is generally satisfied in Bitter pattern experiments on flat 
platelets; see for example, figure (1) in reference [|J. The flux density is determined by, 

B a cos(j) 

B = ~^r- (3) 

In the presence of fixed external currents and temperature the equilibrium flux line lattice 
is obtained by minimising the Gibbs free energy, 

G(H) = l/-5* (4) 

where U contains the interaction between the flux lines and their self-energies, and we have 
defined G in CGS units. 

The Gibbs energy is stationary with respect to the normal component of B, and hence 
the correct thermodynamic potential to be minimised is, 



G = U-^^ (5) 

An v ' 

where Hx is given by equation (2) and Bx by 

B X = B sinO. (6) 

As shown by Kogana, the general expression for the pairwise interaction per unit length, 
U12, between two parallel vortices is 

Un = ^h z (r 12 ) (7) 

where h z (r\2) is the longitudinal component of the magnetic field due to vortex 2 at the 
position of vortex 1 (and vice versa). In a periodic lattice all points are equivalent. Hence, 
the lattice sum is obtained by summing the pairwise interactions between a vortex at the 
origin and another at the lattice position R i5 multiplied by the vortex number density, n. 
Then the interaction energy per unit volume is, 

£/ = ^EMR*)- (8) 

This sum includes the self-energy of the vortices, taken as the limit R; — > 0. 

Using the Poisson sum formula equation (8) can also be expressed as the reciprocal space 
sum, 

tf = ^£*.(Gi)> (9) 

S7r G, 



*j 



where h z (G) is the Fourier transform of h z (H) and Gj are the set of reciprocal lattice 
vectors. 

The prescription for obtaining the equilibrium flux line lattice is now as follows. First U 
is minimised with respect to the vortex positions for a given density and orientation of the 
flux lattice with respect to the c axis. Finally, U is substituted into the Gibbs free energy 
which is minimised as a function of the orientation of the flux lattice. Before U can be 
calculated, however, the longitudinal component of the field distribution must be known. 
This is derived from the London theory in the next section. 



III. ANISOTROPIC LONDON THEORY 

The high temperature superconductors may be viewed as a stack of superconducting 
layers coupled via Josephson tunnelling. This is conveniently described by the Lawrence- 
Doniach modelta The Lawrence-Doniach model introduces four length scales: the in-plane 
coherence length (£ ab ) and penetration depth (X a b), the distance between planes (s) and 



the Josephson length scale, Xj = £ a by 2/p, where p is the dimensionless Josephson coupling 
constant. Xj has the physical significance that for distances from the vortex core of less 
than Xj the phase differences between planes are large. Conversely, for distances from the 
core of greater than Xj the phase differences between planes are small. In the latter case the 
Josephson currents are small and in the extreme type II limit the Lawrence-Doniach model 
becomes equivalent to the three dimensional London theory with uniaxial anisotropy. 

In YBa2Cus07_5 Xj ^60A, which for most typical field strengths is much smaller than 
the inter-vortex spacing. Hence, the vortex-vortex interactions are accurately described 
by London theory. Within the London approximation the basic equation describing the 
magnetic field distribution for an isolated vortex is 

h+(Vx A- Vxh) = <j) 5 (r)z, (10) 

where A = mc 2 /4im s e 2 , and m is the effective mass tensor. For uniaxial anisotropy A has 
two degenerate eigenvalues: A a associated with screening currents flowing in the plane, and 
A c associated with screening currents flowing along the c-axis. The anisotropy ratio, 

7 = VA~ = r (11) 

V A a A a 

is defined by Xj/s. Fourier transforming and inverting equation (10) we obtainO 

h(k) = (TTAS 3 ) 1 Z " I+AJP-MA.-A.WJ *" (12) 

where Q = k x c. 

The field component parallel to the vortex is 
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»,« = Mi±|£l (13) 

where 5 = 1 + X 2 zz k 2 x + X 2 c k 2 , C = 1 + A 2 *; 2 , A; 2 = A; 2 + A; 2 and A 2 Z = A 2 sin 2 # + A 2 cos 2 #. 
Equation (13) is conveniently expressed astil 

h z (k) = k(k) - h 2 (k) = |(| - %^) • (14) 

The expression for h z (r) is obtained from the inverse Fourier transform of h z (k) to produce 

h z {r) = h 1 {r)-h 2 {r), (15) 



where 



h(r) = ^^-Ko(Po) 



r 2 y2 

h — 

A 2 A 2 



Po = ^ + !. (16) 



and 



with 



'^IKVK*^^ (i7 » 



+ I*TT> ( 18 ) 



a 2 (w) 6 2 (m) 
a 2 ( M ) = A 2 z -(AL-A> (19) 



and 



6 2 ( M ) = A 2 -(A 2 -A>. (20) 

fCo and i^i are the zeroth and first order Bessel functions, respectively. 

In isotropic superconductors the magnetic field due to a vortex is parallel to the vortex 
core and positive. Hence, the Lorentz forces between a pair of vortices is centrally directed 
and repulsive. In anisotropic superconductors, however, there are both longitudinal and 



transverse components to the vortex field. Furthermore, the longitudinal component be- 
comes negative with minima situated at a distance x m i n either side of the core in the x 
direction (i.e. in the plane enclosing the c axis and B). The Lorentz force between a pair 
of vortices therefore becomes attractive for certain relative orientations. It is this attrac- 
tive interaction which has led to the suggestion of the vortex-chain regime for flux densities 
where the average vortex spacing becomes comparable to x m i n . In figure 2 we plot x m i n as 
a function of the orientation of the vortex with respect to the c axis for anisotropy ratios of 
5. 

IV. THE LATTICE SUMMATION 

The interaction energy of a lattice of vortices is given in equations (8) and (9). In this 
section we discuss the most efficient way of performing these summations. 

A. Low Field Limit 

In the low field limit, defined by L ^> A, where L is the average inter- vortex spacing, the 
real space sum converges rapidly. The area of the unit cell, s, and hence n, is defined by, 

n B 2 V ; 

The effect of the uniaxial anisotropy is to both lift the orientational degeneracy of the hexag- 
onal lattice and to cause a uniaxial distortion of the equilateral triangles. A deformation of 
(3 X in the x direction and by j3 y in the y direction will produce a family of isosceles triangles. 
For a fixed flux density s is constant, so that (3 x (3 y = 1, and hence we denote (3 X = 1/ ' (3 y = (3. 
There are two orientations of the lattice which are compatible with the uniaxial symmetry. 
However, as shown in ref. [ |3| the orientation of the lattice which minimises the energy for 
high flux densities is that illustrated in figure 3. This is also the orientation which minimises 
the energy in the low density limit. 

A general lattice vector corresponding to the unit cell shown in figure 3 is 
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R-mn = ma.1 + na 2 , 
where ai = L/3x, 

" = L {i A+ W^ (22) 

and m and n are integers. The summations of h±(r) and h2(r) are performed over ellipses 
defined by the contours of ft-i(r) and ^2( r ) 5 respectively. 

In the London approximation the vortex self-energy is divergent due to the inability of the 
theory to describe the core adequately. The theory assumes a constant local magnetic field 
across the core, its value being given at r = £. However, in anisotropic superconductors the 
vortex core is not circular, but elliptical for fields oriented away from c. Hence, an elliptical 
cut-off is required to remove the divergence at small r, namely £ x = £ zz and £ y = £ a , where 
£zz = £a cos2 + £c s i n2 ^- (£a and £ c are the in-plane and out-off-plane coherence lengths) . 
The self-energy is therefore given by 

U se i f - energy = l^x^Ko [-) (23) 

where k = \ ab /£ c = \ c /£ a - 

B. High Field Limit 



For high magnetic fields, L <C A, the real space sums converge slowly. The summation 
of equation (8) can be more easily evaluated in reciprocal space to obtain the lattice energy 

tf = f^£k(G M ), (24) 

«vr pq 

where h z (G pq ) = hi(G M ) - h 2 {G pq ) and 

hl{Gpq) = Al(TTAlMTW (25) 

h(G pq ) = <Po i^A) f 7 rs- (26) 



The coefficients a(u) and b{u) are given in equations (19) and (20). G pq describes the set of 
reciprocal lattice vectors, 

G pq =pbi + qb 2 , 

where bi = G \ — -x y , 

\2(3 2 / 

and b 2 = G /3y, (27) 

with G = 2ttL/s, and p and q are integers. 

The sum of h2(G pq ) converges rapidly. However, the sum of hi(G pq ) is logrithmically 
diverging because of the 1/G 2 behaviour. Physically this corresponds to the self-energy of 
the vortex core. There are various ways to deal with this divergence. One method" is to 
use a Gaussian cut-off in the reciprocal space sum. Another is to use the Ewald summation 
method, originally used by Fetter^ for type II superconductors. The advantage of the latter 
method is that it explicitly removes the formally diverging self energy term, so we adopt it 
in this paper. 

The details of the extension of Fetter's method to anisotropic superconductors are shown 
in the appendix. Here we simply quote the result for the total energy per unit volume, 
excluding the self-energy, as 






2fi2 



exp (— crG 



G 2 G 2 (l + G 2 ) 



1 V E (^ 



Bn<f> fXl - A 2 ^ ^ /-i du (2g) 



8n \ X l J g Jo (l + o?G 2 x + VGl] 
where a 2 = l/4vm, G 2 = A 2 2 G 2 +A 2 G 2 , R 2 = (R rc /A^) 2 + (R y /A c ) 2 and E t is the exponential 
integral. 

The vortex self-energy (obtained in the limit R,; — ► 0) is not included in this summation, 
but can be evaluated as before by employing an elliptical cut-off for the vortex core, namely 
x — > £ x = t; zz and y — > £ y — £ c . Hence, with the substitution Rj = 1/n, the self-energy term 

becomes 

B(f) X zz frm\ 

self-energy 32 ^ X 2\^ { K 2 ) ^ 
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V. RESULTS AND DISCUSSION 

We now turn to a discussion of our predictions for the lattice parameters, and their 
comparsion with experiment. Figure 4a shows our predictions and the experimental values 
from [ I] of the intra-chain (D) and inter-chain (C) distances when projected onto the a-b 
plane, versus the normal component of the applied field. The applied field is at an angle of 40° 
to the c axis. We have taken the anisotropy parameter 7 as 5, the in-plane penetration depth 
A a as 1 413A, and the in-plane coherence length £ a as 16A. For this choice of parameters 
there is very good agreement with experiment. The results for the same parameters are 
shown in figure 4b for an applied field at 70° to the c axis, where the fit is less good. A 
possible explanation for this might be due to a misalignment of the crystal. The inset of 
figure 4b shows the flux line lattice orientation with respect to the c-axis as a function of the 
applied magnetic field. For low magnetic fields the flux lattice is oriented nearly parallel to 
the c axis, slowly aligning with the applied field direction with increasing magnetic intensity. 

To show that the vortex-chain structure has indeed been observed in reference [ [l[ we 
must take into account the effects of the rotation of the flux lattice as a function of applied 
field. This has two consequences: first, the position of the minimum in the vortex- vortex 
interaction, x m i n , will change, as it is angular dependent. This angular dependence is shown 
in figure 2. Second, the magnetic flux density in the plane normal to the flux lattice is 
not linearly related to the applied field, but given by equation (3). We therefore scale the 
results in the following way. The intra-chain distance in the plane normal to the flux lattice, 
d = Dcos9, is divided by x min and plotted against the inverse flux density, 1/B. In the 
vortex-chain regime this quantity should be a constant. To preserve areas the inter-chain 
distance is multiplied by x min , and should be inversely proportional to the flux density. 
Figures 5a and 5b shows the log- log plot of d/x m i n and c x x m i n against 1/B for the applied 
field orientations of 40° and 70°, respectively. We have also plotted the experimental points 
scaled in the same way. Evidently, for low fields, both the theoretical and experimental 
results are in close agreement with the vortex-chain predictions, with d/x m i n roughly unity 
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and almost independent of B. 

At higher fields, scaling of the unit cell parameters vary with the magnetic flux density 
as ~ l/B 1 / 2 , as expected. The anisotropic distortion of the hexagonal lattice agrees with 
the perturbation expansion of [ |3| with, 

^( sin2 * + ; w y 4 . po) 

At high fields the flux lattice is parallel to the applied field, so 9 = <fi. 

As a final comparison to experiment we calculate the value of the intra-chain distance 
on the surface of the sample as a function of the orientation of the applied field for a fixed 
normal component of 12 Oe. This is shown in figure 6 with the experimental values of [ [1| 
for comparison. Again, there is reasonable agreement for this choice of parameters. 

VI. CONCLUSIONS 

In summary, we have calculated the equilibrium flux line lattice in planar crystals of 
YBa 2 Cu307_,5 using three dimensional anisotropic London theory. By taking into account 
demagnetisation effects, which cause the flux line lattice to orient away from the applied field 
towards the crystalline axis, we are able to show that the low field Bitter pattern experiments 
of Gammel et alB demonstrate the existence of the "vortex-chain" state. In this state the 
intra-chain distance is independent of flux density, whereas the inter-chain distance scales 
as the inverse flux density, 1/B. We also demonstrate that as the field strength is increased 
there is a smooth cross-over to the distorted hexagonal lattice in which the inter-vortex 
spacings scale uniformly as 1/B 1 / 2 . 

The agreements between theory and experiment are found for an anisotropy ratio 7 of 
5 and the in-plane penetration depth A a of 1 413A. This value of the penetration depth is 
the zero temperature penetration depth, and not the value of the penetration depth at the 
irreversibility line quoted for twinned samples. Since the samples used in ref [ Q are relatively 
twin-free, consisting of untwinned regions of at least 100/zm square, we deduce that locally 
the lattice is not frozen but has assumed its zero temperature equilibrium configurations. 
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APPENDIX A: 



In this appendix we use the Ewald summation method to derive the interaction energy 
of a lattice of vortices in the high density limit. 
The reciprocal space sum 

1 



s=J2 



allG 



1 + G 2 



(Al) 



diverges logarithmically for large G, where G 2 = \ 2 ZZ G 2 X + A 2 G 2 . 
(Al) may be written as 



g^o G2 ffoGMlfG 2 



(A2) 



where the third term converges, leaving the second term to consider. Using the identity 



&= 2 J £exp(-£ 2 G 2 )^, 
and splitting the integral into two parts, we obtain 



(A3) 



E — = E 

G^O G^O 

= E 

G^O 



i exp -e& )d£+ £ exp -£ 2 G 2 d£ 



-2/~i2 



2y o eexp(-e 2 G^de 



2/^2 



exp (— a G 
G 2 



(A4) 



where the choice of a = (Ann)^ 1 ^ 2 maximises the rate of convergence and n is the real space 
density of lattice points. The second term of (A4) now has a Gaussian cut-off and converges 
easily. However, the logarithmic divergence has to be extracted from the first term. Writing 
the first term of (A4) as a sum over all G, 



T = E 2 ^exp(-£ 2 G 2 )^-2^^ 



allG 

and using the Poisson sum formula, n^Z G F{G) = Z)rF(R), we obtain 

1 

allRi 



T 



4nnX zz X c ^ 2 L^\-—2)- 



5l\ ^ 



a 



(A5) 
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where R 2 = (R x /\ zz ) 2 + (Rj,/A c ) 2 . 
Finally, collecting terms, 



2/^2 



exp ( — a G 



G 2 1 + G 2 






— a 



where the exponential integral is 



(A6) 



00 exp(— i) 



El ( X ) = r e -^>dt. 

Jx t 



The logarithmically diverging contribution comes from the term Rj — > in (A6). 
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Figure Captions 

Figure 1. The slab geometry used in the decoration experiments, with the c-axis per- 
pendicular to the surface of the slab. The (X, Y, Z) axes define the crystallographic axes 
with Z being parallel to c and (X — Y) lying in the (a — b) plane. The external field, H a , is 
applied at an angle to c, producing a vortex lattice of magnetic flux density, B, inclined 
at an angle 9 to c. This vortex lattice is described by a general coordinate axes (x, y, z) by 
a rotation of 9 about the Y-axis such that c, x and z are coplanar and B is parallel to z. 

Figure 2. The position of minimum of the inter- vortex potential, x m i n , in units of A = 
VAaAc versus the orientation of the vortex for 7 = 5. 

Figure 3. (a) The isotropic flux line lattice is described by equilateral triangles of side 
L. (b) A deformation of (3 in the x direction and by 1/(3 in the y direction will produce a 
family of isosceles triangles that describe the anisotropic flux line lattice. 

Figure 4(a). Intra-chain distance (D): theory, solid line; experiment, open diamonds, and 
inter-chain distance (C): theory, dashed line; experiment, filled diamonds, projected onto the 
a-b plane. Applied field orientation of 40°, in-plane penetration depth A a of 1 413A, in-plane 
coherence length £ a of 16A and anisotropy ratio 7 of 5. The inset shows the vortex chain. 

Figure 4(b). The same as fig 4(a) with an applied field orientation of 70°. The inset 
shows the orientation of the flux line lattice as a function of the applied field: solid line - 
0=40° and dotted line - 0=70°. 

Figure 5(a). Log-log plot of intra-chain distance divided by x m i n : theory, solid line; 
experiment, open diamonds and inter-chain distance multiplied by x m i n : theory, dashed 
line; experiment, filled diamonds, in the plane normal to the flux lattice. Applied field at 
40°. 

Figure 5(b). The same as fig 5(a) with the applied field at 70°. 

Figure 6. The intra-chain distance on the surface of the sample as a function of on the 
applied field orientation for a fixed normal component of 12 Oe, open squares - experiment. 
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